o 
o 

(N 

o 
O 

0^ 



BELLE-CONF-0916 
draft v.4.0 

Measurement of the form factors of the decay B'^ — ^ D*^i'^Ui and 
determination of the CKM matrix element \Vcb\ 

I. Adachi/° H. Aihara,^^ K. Arinstein/' T. Aso,''^ V. Aulchenko/'^^ T. Aushev,26.i9 

T. Aziz,^^ S. Bahinipati,^ A. M. Bakich,^^ y Balagura/^ Y. Ban,^^ E. Barberio,^^ 
A. Bay,26 I. Bedny/-^^ K. Belous/^ V. Bhardwaj,^^ B. Bhuyan/^ M. Bischofberger,^^ 
S. Blyth,35 A. Bondar/'^i A. Bozek,^^ M. Bracko,28'20 j. Brodzicka,^^ T. E. Browder,^ 
M.-C. Chang,^ P. Chang,^^ Y.-W. Chang,^^ Y. Chao,='^ A. Chen,^^ K.-F. Chen,^^ 
P.-Y. Chen,36 B. G. Cheon,^ C.-C. Chiang,^^ R. Chistov/^ I._S. Cho,^^ S.-K. Choi/ 
Y. Choi,52 J Cmkovic/2 J. Dalseno,29'55 M. Danilov/^ A. Das,^^ M. Dash,^^ 

A. Drutskoy,^ W. Dungel/^ S. Eidelman/'^^ D. Epifanov/'^^ M. Feindt,^^ H. Fujii/° 
M. Fujikawa, N. Gabyshev/' "'^ A. Garmash/'^^ G. Gokhroo,^"^ P. Goldenzweig,^ 

B. Golob,27.20 M. Grosse Perdekamp/2.47 g 49 ^ 23 j Haba/° B.-Y. Han,^^ 

K. Hara,^^ T. Hara/° Y. Hasegawa,^^ N. C. Hastings, K. Hayasaka,^^ H. Hayashii,^^ 
M. Hazumi,i° D. Heffeman,^^ rj Higuchi,i° Y. Horii,^^ Y. Hoshi,^^ K. Hoshina,'^^ 
W.-S. Hou,3^ Y. B. Hsiung,36 H. J. Hyun,^^ Y. Igarashi,i° T. lijima,^^ K. Inami,^! 
. Ishikawa,^^ H. Ishino,^°'|!| K. Itoh,^^ R. Itoh,^° M. Iwabuchi,*^ M. Iwasaki,^^ Y. Iwasaki, 
\ Jinno,3i M. Jones,^ N. J. Joshi,^^ T. Julius,^^ D. H. Kah,^^ H. Kakuno,^^ J. H. Kang,^ 



H 
R 

M, 



0^ . 

O ■ R 



X 



i\. isiiiKawa, n. ismiio, iv. ixoii, sx. ixoii, m. iwaDucm, m. iwasaKi, i. iwasam, 
T. Jinno,3i M. Jones,^ N. J. Joshi,^^ T. Julius,^^ D. H. Kah,^^ H. Kakuno,^^ J. H. Kang,^^ 
P. Kapusta,^'' S. U. Kataoka,^^ N. Katayama,^" H. Kawai,^ T. Kawasaki, A. Kibayashi,^° 
H. Kichimi,i° C. Kiesling,^^ H. J. Kim,^^ H. O. Kim,^^ J. H. Kim,^^ g Kim,^^ 
Y. 1. Kim,25 Y. J. Kim,6 K. Kinoshita,^^ B. R. Ko,^^ S. Korpar,28'20 M. Kreps,^^ 
P. Krizan,27'20 p Krokovny,io T. Kuhr,22 R. Kumar,^^ T. Kumita,'^^ E. Kurihara,^ 
E. Kuroda,^^ Y. Kuroki,^^ Kusaka,^^ A. Kuzmin,i'^^ Y.-J. Kwon,^^ S.-H. Kyeong,^^ 
J. S. Lange,^ G. Leder,^^ M. J. Lee,^° S. E. Lee,^° S.-H. Lee,^^ J. Li,^ A. Limosani,^^ 
S.-W. Lin,36 C. Liu,^9 D. Liventsev,^^ R. Louvot,^^ J. MacNaughton,i° F. Mandl,^^ 
Marlow,^*^ A. Matyja,^^ S. McOnie,^^ T. Medvedeva,^^ Y. Mikami,^^ K. Miyabayashi,^^ 
Miyake,^^^ H. Miyata,^^ Y. Miyazaki,^! R. Mizuk,^^ A. Moll,29'55 t. Mori,^! T. Miiller,22 
Mussa,^® T. Nagamine,^^ Y. Nagasaka,"^"*^ Y. Nakahama,^^ I. Nakamura,"^" E. Nakano,"^^ 
Nakao,^° H. Nakayama,^^ H. Nakazawa,^^ Z. Natkaniec,^^ K. Neichi,^^ S. Neubauer,^^ 
S. Nishida,i° K. Nishimura,^ O. Nitoh,'52 g_ Noguchi,^^ T. Nozaki,^^ A. Ogawa,^^ 
3. Ogawa,56 T. Ohshima,^! S. Okuno,^! S. L. 01sen,5° W. Ostrowicz,^^ H. Ozaki,^^ 
P. Pakhlov,i9 G. Pakhlova,i9 H. Palka,^^ C. W. Park,52 g Park,25 H. K. Park,25 
K. S. Park,^2 l g_ Peak,^^ Pernicka,^^ R. Pestotnik,^^ M. Peters,^ L. E. Piilonen,^^ 
A. Poluektov,^''^^ K. Prothmann,29' b. Riesert,^^ M. Rozanska,^^ H. Sahoo,^ 
K. Sakai,^^ Y. Sakai,^° N. Sasao,^^ O. Schneider, 2^ P. Schonmeier,^*^ J. Schiimann,^^ 
C. Schwanda,!^ A. J. Schwartz,^ R. Seidl,^^ A. Sekiya,^^ K. Senyo,^! M. E. Sevior,^^ 
L. Shang,i4 M. Shapkin,!^ V. Shebahn,!'^^ C. P. Shen,^ H. Shibuya,^^ g_ ghiizuka,^! 
S. Shinomiya,^^ j _q shiu,^^ B. Shwartz,!'^^ F. Simon,29'55 ^ Singh,^^ R. Sinha,^^ 
A. Sokolov,i6 E. Solovieva,i9 S. Stanic,^° M. Staric,^^ J. Stypula,^^ A. Sugiyama,^^ 
K. Sumisawa,^° T. Sumiyoshi,^^ S. Suzuki,^^ S. Y. Suzuki, Y. Suzuki, F. Takasaki,^'' 
N. Tamura,39 K. Tanabe,^^ M. Tanaka,^° N. Taniguchi,^° G. N. Taylor,30 Y. Teramoto,^^ 
I. Tikhomirov,^^ K. Trabelsi,^^ Y. F. Tse,^° T. Tsuboyama,^° K. Tsunada,^^ Y. Uchida,^ 

S. Uehara,^° Y. Ueki,^! K. Ueno,^^ T. Uglov,^^ y. Unno,^ S. Uno,i° P. Urquijo,^^ 
Y. Ushiroda,^° Y. Usov,^'^^ G. Varner,^ K. E. Varvell,^^ K. Vervink,^^ A. Vinokurova,^''^^ 



Typeset by REVI^ 



C. C. Wang,36 C. H. Wang,^^ J. Wang,^^ M.-Z. Wang,^^ P. Wang/^ X. L. Wang/^ 
M. Watanabe,39 Y. Watanabe,^! R. Wedd,^" J.-T. Wci,-'^^ J. Wicht/o L. Widhalm/^ 
J. Wiechczynski,^'' E. Won,^^ B. D. Yabsley,^^ H. Yamamoto;^^ Y. Yamashita,^^ 
M. Yamauchi/° C. Z. Yuan/^ Y. Yusa,^^ C. C. Zhang/^ L. M. Zhang,^^ Z. P. Zhang,^^ 
V. Zhilich/'^i V. Zhulanov/'^^ T. Zivko,^^ A. Zupanc,^^ N. Zwahlen,^^ and O. Zyukova^'^^ 

(The Belle Collaboration) 

(The Belle Collaboration) 

^Budker Institute of Nuclear Physics, Novosibirsk 

^Chiha University, Chiha 
^University of Cincinnati, Cincinnati, Ohio 45221 
^Department of Physics, Fu Jen Catholic University, Taipei 
^ Justus- Liebig-Universitat Giefien, Giefien 
^The Graduate University for Advanced Studies, Hayama 
'^Gyeongsang National University, Chinju 
^Hanyang University, Seoul 
^University of Hawaii, Honolulu, Hawaii 96822 
^^High Energy Accelerator Research Organization (KEK), Tsukuba 

^^Hiroshima Institute of Technology, Hiroshima 
University of Illinois at Urbana- Champaign, Urbana, Illinois 61801 
^^India Institute of Technology Guwahati, Guwahati 
^'^ Institute of High Energy Physics, 
Chinese Academy of Sciences, Beijing 
^^Institute of High Energy Physics, Vienna 
Institute of High Energy Physics, Protvino 
^'^ Institute of Mathematical Sciences, Chennai 
^^INFN - Sezione di Torino, Torino 
^^Institute for Theoretical and Experimental Physics, Moscow 
^^J. Stefan Institute, Ljubljana 
^^Kanagawa University, Yokohama 
^^Institut fiir Experimentelle Kernphysik, Universitdt Karlsruhe, Karlsruhe 

^^Korea University, Seoul 
^'^Kyoto University, Kyoto 
Kyungpook National University, Taegu 
^^Ecole Polytechnique Federale de Lausanne (EPFL), Lausanne 
^'^ Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana 

University of Maribor, Maribor 
Max-Planck-Institut fiir Physik, Miinchen 
University of Melbourne, School of Physics, Victoria 3010 
^^Nagoya University, Nagoya 
^^Nara University of Education, Nara 

^^Nara Women's University, Nara 
^^National Central University, Chung-li 
^^National United University, Miao Li 
^^Department of Physics, National Taiwan University, Taipei 
^^H. Niewodniczanski Institute of Nuclear Physics, Krakow 



2 



^^Nippon Dental University, Niigata 
^^Niigata University, Niigata 
University of Nova Gorica, Nova Gorica 
^^Novosibirsk State University, Novosibirsk 
Osaka City University, Osaka 
^^Osaka University, Osaka 
'^^Panjab University, Chandigarh 
^^Peking University, Beijing 
Princeton University, Princeton, New Jersey 08544 
^''RIKEN BNL Research Center, Upton, New York 11973 
'^^Saga University, Saga 
University of Science and Technology of China, Hefei 
^'^Seoul National University, Seoul 

^^Shinshu University, Nagano 
^^Sungkyunkwan University, Suwon 
School of Physics, University of Sydney, NSW 2006 
Tata Institute of Fundamental Research, Mumbai 
Excellence Cluster Universe, Technische Universitdt Miinchen, Garching 
^'^Toho University, Funahashi 
^'^Tohoku Cakuin University, Tagajo 
^^Tohoku University, Sendai 
Department of Physics, University of Tokyo, Tokyo 
Tokyo Institute of Technology, Tokyo 
Tokyo Metropolitan University, Tokyo 
Tokyo University of Agriculture and Technology, Tokyo 
^^Toyama National College of Maritime Technology, Toyama 
^'^IPNAS, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061 

^^Yonsei University, Seoul 

Abstract 

We present a measurement of the Cabibbo-Kobayashi-Maskawa matrix element \Vch\ times the 
form factor normalization J-{1) using the decay — > D*^£'^v^. This measurement is performed 
together with a determination of the form factor parameters p^, -Ri(l) and -R2(l) which fully 
characterize this decay in the framework of Heavy Quark Effective Theory. This analysis is based 
on a data sample equivalent to 140 fb~^ of Belle data collected near the T(4S) resonance. 

The preliminary results, based on about 27,000 reconstructed B^ D*^i~^V(i decays, are = 
1.376 ±0.074 ±0.056, i?i(l) = 1.620±0.091±0.092, i?2(l) = 0.805 ±0.064 ±0.036 and J^(l)|V;b| = 
35.0±0.4±2.2. We find the B+ D*°£+Ui branching fraction to be (4.84±0.04±0.56)%. For all 
numbers quoted here, the first error is the statistical and the second is the systematic uncertainty. 

A direct, model-independent determination of the form factor shapes has also been carried out 
and shows good agreement with the HQET based form factor parametrization by Caprini et ai. 
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FIG. 1: Quark-level Feynman diagram for the decay D*^t^u^. 

I. INTRODUCTION 

The magnitude of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element \Vcb\ can 
be determined from exclusive decays B D*ii' [l]. In the framework of Heavy Quark 
Effective Theory (HQET), this decay is described by three form factor parameters, p^, 
Ri{l) and -R2(l) [21, y] along with a form factor normalization J-'{1) derived from lattice 
QCD Many measurements of the mode are available which are unfortunately barely 
consistent [sj, suggesting a hidden systematic uncertainty. 

To shed light on this inconsistency, we now investigate the related decay B~^ D*^t^v^ 
which carries different experimental systematics. This is interestirig also because only few 
measurements of the B^ mode are available in the literature jl, 0, HI • 

This paper is organized as follows: after introducing the theoretical framework for the 
study of this decay, the experimental procedure is presented in detail. This is followed by 
a discussion of our results and the systematic uncertainties. Finally, a form factor parame- 
terization independent measurement of the form factor shapes is described and results are 
given. 




II. THEORETICAL FRAMEWORK 
A. Kinematic variables 



The decay B^ D*^i~^U£ [9| proceeds chiefly through the tree-level transition shown in 
Fig. [H Its kinematics can be fully characterized by four variables: 
The flrst one is w, deflned by 

Pb ■ Pd*o ml + m^.o - s 
w = = , (Ij 

mBfTiDto 2mBm£)*o 

where ttlb and m/j.o are the masses of the B~^ and the D*^ mesons (5.2792 GeV/c^ and 
2.007 GeV/c^, respectively [iC?]), Pb and Pd*o are their four-momenta, and = (P^ + P,^)^. 
In the B rest frame, approximately equal to the T(45') center-of-mass (cm.) frame, the 
expression for w reduces to the Lorentz boost Jd*o = _E'£)*o/m^*o. The ranges of w and 
are restricted by the kinematics of the decay, with = corresponding to 

W„,ax = ^ — ~ 1.505 , 2 

and Wmin = 1 to 

^Lx = {rriB - niD^o)' ^ 10.71 GeV'/c' . (3) 
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FIG. 2: Definition of the angles Oi, By and x for the decay B+ D*^i+ue, ^ D^tt'^^. 

The point w = 1 is also refered to as zero recoil. 

The remaining three variables are the angles shown in Fig. [2) 

• 6i, the angle between the direction of the lepton in the virtual W rest frame and the 
direction of the W in the B rest frame; 

• 6y, the angle between the direction of the D meson in the D*^ rest frame and the 
direction of the D*^ meson in the B rest frame; 

• X) the angle between the plane formed by the D*^ and the plane formed by the W decay. 



B. Four-dimensional decay distribution 

The Lorentz structure of the B^ — > D*^i~^i'£ decay amplitude can be expressed in terms 
of three helicity amplitudes {H^, and Hq), which correspond to the three polarization 
states of the D*'^, two transverse and one longitudinal. For low- mass leptons (electrons and 
muons), these amplitudes are expressed in terms of the three functions /i^j(w), Ri{w), and 




Hi{w) = rriB 



R*{l-r^){w + l) 



=-hAAw)Hi{w) , 



(4) 



2 Vl - 2wr + r2 
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where 



Vl-2w + r2 (l ± ,J^,Ri{w] 



1 + 



1 — r 
{w - - R2{w)) 



(5) 
(6) 



with R* = {2yJmBmj;)*o) / {niB + m£)*o) and r = m£,*o/mB- The functions -Ri(w) and R2{w) 
are defined in terms of the axial and vector form factors as, 



R,iw) 



V{w) 



w+1 

Ri{w) 2 
w + 1 

By convention, the function (w) is defined as 

1 2 



Ai{w) . 



R*w + 1 



AH 



(7) 
(8) 

(9) 



For w 1, the axial form factor Ai{w) dominates, and in the limit of infinite b- and c-quark 



masses, a single form factor describes the decay, the so-called Isgur-Wise function Il2 
The fully differential decay rate in terms of the three helicity amplitudes is 



d^r(5+ ^ D*H+ue 



■Vw^ - 1(1 - 2wr + r^)Gl\V, 



d^i;d(cose^)d(cos^v')dx 8{Att)^ 
X {{1 ~ cosOe)^ sin^ OvHliw) + {1 + cos 6 sinHy H'^_{w) 

+ 4 sin^ Be cos^ OyH^iw) - 2 sin^ 6^ sin^ By cos 2xH+{w)H_{w) 
— 4sin6'^(l — cos^^^) sin^^y cos6'y cos xH+{w)Hq{w) 

+ 4sin^£(l + cos^^£) sin6'y cos6'v' cos xH-{w)Ho{w)^ , 



cb\ 



(10) 



with Gf = (1.16637 ± 0.00001) x 10"^ GeV-^ [lO|. By integrating this decay rate over all 
but one of the four variables, w, cos 9i, cosOy, or x, we obtain the four one-dimensional 
decay distributions from which we will extract the form factors. The differential decay rate 
as a function of w is 



dr 



dw 487r3 



q2 ^ 

^ m\,,o{mB - niD*o) Q{w)J^'^{w)\Vcb\^ 



(11) 



where 



T^{w)g{w) = /ii^(u;)V^^7^(w + 1)2 J 2 



1 — 2wr + r 



X 



l + Ri(w) 



,w — 1 
w + 1 



l + il-R2{w)) 



'1 -r 



w 
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and Q{w) is a known phase space factor, 



In the infinite quark-mass hmit, the heavy quark symmetry (HQS) predicts = 
1. Corrections to this hmit have been calculated in lattice QCD. The most recent result 
obtained in unquenched lattice QCD reads J^(l) = 0.921 ± 0.013 ± 0.020 

C. Form factor parameterization 

The heavy quark effective theory (HQET) allows a parameterization of these form-factors 
to be obtained. Perfect heavy quark symmetry implies that Ri{w) = R2{w) = 1, i.e., 
the form factors A2 and V are identical for all values of w and differ from Ai only by a 
simple kinematic factor. Corrections to this approximation have been calculated in powers 
of AqcD/'^b and the strong coupling constant as- Various parameterizations in powers of 
{w — 1) have been proposed. Among the different predictions relating the coefficients of the 
higher order terms to the linear term, we adopt the following expressions derived by Caprini, 
Lellouch and Neubert [3|, 

Ha, (w) = hA,{l) [1 - 8p^z + (53p^ - 15)z^ - {231p^ - 91)/] , (12) 
Ri{w) = Ri{l) - 0.12(w - 1) + 0.05(w - 1)^ , (13) 
R2(w) = i?2(l) + 0.11{w - 1) - 0.06(w - If , (14) 

where z = {^/w + 1 — \/2)/{-\/w + 1 + \/2). The three parameters p^, -Ri(l), and -R2(l), 
cannot be calculated; they must be extracted from data. 

III. EXPERIMENTAL PROCEDURE 
A. Data sample and event selection 

The data used in this analysis were taken with the Belle detector [13] at the KEKB 
asymmetric energy e^e~ collider [3]. Belle is a large-solid-angle magnetic spectrometer 
that consists of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), 
an array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrangement of time- 
of-fiight scintillation counters (TOF), and an electromagnetic calorimeter (ECL) comprised 
of CsI(Tl) crystals located inside a super-conducting solenoid coil that provides a 1.5 T 
magnetic field. An iron flux-return located outside of the coil is instrumented to detect 
mesons and to identify muons (KLM). In the original setup the SVD was composed of three 
layers of double sided silicon strip modules, in 2003 it was upgraded and a fourth layer was 
added. 

The data sample consists of 140 fb^^ taken at the T(4S') resonance, or 152 x 10^ 
BB events. The entire data sample has been recorded using the three-layer SVD setup. 
Another 15 fb~^ taken at 60 MeV below the resonance are used to estimate the non-BB 
(continuum) background. The off-resonance data is scaled by the integrated on- to off- 
resonance luminosity ratio corrected for the 1/s dependence of the qq cross-section. This 
sample is identical to the one used for the B^ D*^£~^Ui analysis |l5| . 



1 + 4- 



w 



1 — 2wr + r 



w + 1 1 
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Monte Carlo (MC) samples equivalent to about three times the integrated luminosit 
are used in this analysis. MC simulated events are generated with the evtgen program 
and full detector simulation based on GEANT [17|] is applied. QED bremsstrahlung in 
B — > Xiiy decays is added using the PHOTOS package [l8|. 

Hadronic events are selected based on the charged track multiplicity and the visible 
energy in the calorimeter. The selection is described in detail elsewhere [l9|. We also apply 
a moderate cut on the ratio of the second to the zeroth Fox- Wolfram moment 13], R2 < 0.4, 
to reject continuum events. 



B. Event reconstruction 

Charged tracks are required to originate from the interaction point by applying the fol- 
lowing selections on the impact parameters in R(j) and z, dr < 2 cm and \dz\ < 4 cm, 
respectively. Additionally, we demand at least one associated hit in the SVD detector. For 
pion and kaon candidates, the Cherenkov light yield from ACC, the time-of-flight informa- 
tion from TOF and dE/dx from CDC are required to be consistent with the respective mass 
hypothesis. 

Neutral D meson candidates are searched for in the decays channels — >■ K~Tt^ and 
— > K'TT^TT'TT^ . We fit the charged tracks to a common vertex and reject the candi- 
date if the x^-probability is below 10^'^. The momenta of the charged tracks are re-evaluated 
at the vertex and the 4-momentum is calculated as their sum. The reconstructed mass 
is required to lie within ±3 standard deviations from the nominal mass [lo|, where one sigma 
(measured from real data) is about 4.5 MeV/c^ (4 MeV/c^) for the one (three) pion mode. 

The D*^ candidate is combined with a slow neutral pion vr^ to form a D*^ candidate. The 
7r° is obtained by combining two photons with > 100 MeV for 9 < 32°, E^ > 150 MeV 
for 6 > 130°, and E^ > 50 MeV for photons detected in the barrel region. The invariant 77 
mass has to lie within 3 standard deviations of the nominal vr*^ mass, with a ~ 5 MeV/c^. 
Additionally, a mass constraint fit is performed and the fit probability has to exceed 10~^. 
Continuum suppression is achieved by requiring a -D*° momentum less than 2.45 GeV/c in 
the cm. frame. 

Finally, the D*^ candidate is combined with a lepton (electron or muon), appropriately 
charged with respect to the kaon. Electron candidates are identified using the ratio of 
the energy detected in the ECL to the track momentum, the ECL shower shape, position 
matching between track and ECL cluster, the energy loss in the CDC and the response of 
the ACC counters. Muons are identified based on their penetration range and transverse 
scattering in the KLM detector. In the momentum region relevant to this analysis, charged 
leptons are identified with an efficiency of about 90% and the probability to misidentify a 



pion as an electron (muon) is 0.25% (1.4%) [2l|, |22]. No SVD hit requirement is made for 



lepton tracks. In the lab frame, the (transverse) momentum of the lepton is required to 
exceed 0.80 GeV/c (0.65 GeV/c) in case of electrons and 0.85 GeV/c (0.75 GeV/c) in case 
of muons. We also apply an upper lepton momentum cut at 2.4 GeV/c in the cm. frame to 
reject continuum. 

In electron events, we attempt bremsstrahlung recovery by searching for photons in a 
cone of 3° around the electron track. If such a photon is found it is merged with the electron 
and the sum of the momenta is assumed to be the lepton momentum. 

Figs. [3] and H] show the invariant mass of the candidates and the 7r° momentum 
distributions, respectively. 
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FIG. 3: Invariant candidate mass distributions in the different sub-samples. All analysis cuts 
(except on the plotted variable) are applied. The estimation of the background contributions is 
described in section IIII CI 

C. Background estimation 

Because we do not reconstruct the other B meson in the event, the B momentum is a 
priori unknown. However, in the cm. frame, one can show that the B direction hes on a 



cone around the {D i)-axis [23 



cose B,D*H = oi^ 11^ 1 ' (15) 



^\Ph\\P 



— *): 



where E'^ is half of the cm. energy and |jr^| is \/ E*^ — m?^. The quantities -Z?^,o^, P/j.o^ and 
mD*oi are calculated from the reconstructed D*^i system. 

This cosine is also a powerful discriminator between signal and background: signal events 
should strictly lie in the interval (—1, 1), although - due to finite detector resolution - about 
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p [GeV/c ] p „, [GeV/c ] 



FIG. 4: Momentum distributions of the slow tt^ (lab frame) in the different sub-samples. All 
analysis cuts are applied. The color scheme is explained in Fig.[3l The estimation of the background 
contributions is described in section [III CI 

5% of the signal is reconstructed outside this interval. The background on the other hand 
does not have this restriction. 

To determine the amount of signal and background in the selected events, we perform a 
two dimensional binned fit to cos 63,0*01 and Am, the invariant mass difference between the 
and the D candidates. We consider backgrounds from the following six sources: 

• continuum: any candidate reconstructed in a non-T(4S') event, 

• fake D^: the candidate has been misreconstructed, 

• combinatoric D*^: the D*^ candidate is misreconstructed, however the candidate 
is identified correctly, 

• fake lepton: the lepton has been misidentified but the D*^ is reconstructed properly, 

• uncorrelated background: the D*^ and the lepton stem from different B mesons, 

• D**: background from B decays with D** D*^n 01 B ^ D*°iTi'^u 
non-resonant, and 

• correlated background: background from other processes in which the Z)*° and the 
lepton stem from the same B meson, e.g, B^ —>■ D*^t^u, r+ yU+z/z/. 
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FIG. 5: Result of the fits to the cosO^ £)*0£ vs. Am distributions in the KTr,e sub-sample. The 
projections in cosO^ £,*oi and Am are shown. The color scheme is explained in Fig. [3l 




FIG. 6: Same of Fig. [5] for the KTT,fi sub-sample. 



These background components are modeled by MC (Z}*°£) events, appropriately selected 
using generator information, except for continuum which is modeled by off-resonance events. 
For muon events, the shape of the fake lepton background is corrected by the ratio of the 
pion fake rate in the experimental data over the same quantity in the MC. The lepton 
identification efficiency is reweighted to compensate for imperfect MC simulation of the 



detector. The actual fit uses the TFractionFitter algorithm [2^ in ROOT j25|]. The fit 
is done separately in the four sub-samples defined by the decay channel and the lepton 
type. The results are shown in Figs. E] to [Hand Table [H 

In all fits, the continuum normalization is fixed to the on- to off-resonance luminosity 
ratio, corrected for the 1/s dependence of the e^e~ qq cross-section. The on- to off- 
resonance luminosity ratio is measured using Bhabha events to an accuracy of 1.5 %. Aside 
from the before mentioned corrections to the lepton identification, the normalization of the 
fake lepton component has also been fixed to the MC expectations, since the fit failed to 
determine this small background contribution reliably. For the uncertainty in this component 
a very conservative value of 30% is assumed. The determination of the D** and signal 
correlated backgrounds works reliably in the two — > K~7i'~^ modes, but fails in the three 
pion channels. In these cases we fix the data to MC ratio to the value found in the former 
channels. We assume the relative uncertainty found in the Kir modes. 
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FIG. 8: Same of Fig. [5] for the K3-K,fi sub-sample. 





Ktt, e 


Ktt, fi 


K37r, e 


i^Syr, fj, 


Raw yield 


13035 


12262 


16989 


16350 


Signal events 


8133 ± 205 


7447 ± 201 


5987 ± 229 


5539 ± 222 


Signal 


(62.39 ± 1.57)% 


(60.73 ± 1.64)% 


(35.24 ± 1.35)% 


(33.88 ± 1.36)% 


Signal correlated 


(1.27 ± 0.31)% 


(1.46 ± 0.32)% 


(1.16 ± 0.26)% 


(1.34 ± 0.31)% 




(0.77 ± 0.98)% 


(0.73 ± 0.98)% 


(0.39 ± 0.50)% 


(0.36 ± 0.47)% 


Uncorrelated 


(4.97 ± 0.54)% 


(4.25 ± 0.45)% 


(3.48 ± 0.41)% 


(3.30 ± 0.38)% 


Fake i 


(0.31 ± 0.10)% 


(1.94 ± 0.59)% 


(0.18 ± 0.06)% 


(0.95 ± 0.29)% 


Combinatoric D*^ 


(24.76 ± 0.51)% 


(24.30 ± 0.48)% 


(16.35 ± 0.69)% 


(15.19 ± 0.67)% 


Fake D° 


(2.91 ± 0.25)% 


(3.12 ± 0.23)% 


(38.53 ± 0.50)% 


(39.45 ± 0.51)% 


Continuum 


(2.63 ± 0.43)% 


(3.46 ± 0.51)% 


(4.68 ± 0.50)% 


(6.14 ± 0.56)% 



TABLE I: The signal and background fractions for selected events within the signal window 
I coseB^D-oel < 1 and 140 MeV < Am < 144.5 MeV. 
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FIG. 9: Reconstruction of the B'^ direction. Refer to the text for details. 

In Table H] those components which have been fixed in the fit are shown with the corre- 
sponding error estimate described above. The error of the continuum component is appro- 
priately scaled by the on- to off-resonance luminosity ratio. 

In the rest of this analysis, we require D*^i candidates to pass the signal window require- 
ment -1 < cos 9B,D'0e < 1 and 140 MeV < Am < 144.5 MeV. 

D. Kinematic variables 

To calculate the four kinematic variables - w, cos 6i, cosOy and x ^ that characterize 
the D*^l^v decay defined in Sect. Ill At we need to determine the B rest frame. The 

B direction is already known to be on a cone around the (Z)*°£)-axis with opening angle 
26B,D*oe in the cm. frame, Eq. [151 For the best determination of the B direction, we first 
estimate the cm. frame B vector by summing the momenta of the remaining particles in 
the event (p^nciusive 01) choose the direction on the cone that minimizes the difference 
to Prelusive, Fig- M 

To obtain p^nciusivo ^^^t apply the following track and cluster selection requirements 
to particles which have not been associated with the {D*^i) system: neutral clusters are 
required to satisfy > 100 MeV for polar angles with respect to the beam direction 
9 < 32°, > 150 MeV for 9 > 130° and E^ > 50 MeV in the barrel region in between; 
the impact parameter cuts for charged tracks are dr < 20 cm and \dz\ < 100 cm for 
Pt < 0.25 GeV/c, dr < 15 cm and \dz\ < 50 cm for px < 0.5 GeV/c, and dr < 10 cm and 
\dz\ < 20 cm for > 0.5 GeV/c. We exclude duplicated tracks by selecting pairs of charged 
tracks with < 275 MeV/c, momentum difference Ap < 100 MeV/c and relative angle 
smaller than 15° or larger than 165°. If such a pair is found, only the track with smaller 
value of (5 dr)^ + {\dz\Y is retained. 

Then, we compute ^inclusive (in the lab. frame) by summing the 3-momenta of the selected 
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particles, 

^inclusive = PHER + PLER -^Pi , (16) 

i 

where the indices HER and LER correspond to the coUiding beams, and transform this 
vector into the cm. frame. Note that we do not make any mass assumption for the charged 
particles. The energy component of the 4-vector Pinciusive is defined by requiring -Einciusivo 
be E^^^^ = v/i/2. 

With the B rest frame reconstructed in this way, the resolutions in the kinematic variables 
are found to be about 0.021, 0.049, 0.057 and 6.46° for w, cos9e, cos9v and x, respectively. 



E. Fit procedure 

We perform a binned fit of the w, cos 9i, cosOy and x distributions over (almost) 
the entire phase space to measure the following quantities: the form factor normalization 
jF(l)|V^b|, Eq- El and the three parameters p^, -Ri(l) and -R2(l) which parameterize the 
form factor in the HQET framework, Eqs. [T2] - [T^ Instead of fitting in four dimensions, we 
fit the one-dimensional projections of w, cos6e, cos 6'y and x to have enough entries in each 
bin of the fit. This introduces bin-to-bin correlations which have to be accounted for. 

The distributions in w, cos9e, cos 6'y and x are divided into ten bins of equal width. The 
kinematically allowed values of w are between 1 and ~ 1.505 but we restrict the w range 
to values between 1 and 1.5. In each sub-sample, there are thus 40 bins to be used in 
the fit. In the following, we label these bins with a common index i, i = 1, . . . ,40. The 
bins i = 1, . . . , 10 correspond to the bins of the w distribution, i = 11, . . . , 20 to cos 9i, 
z = 21, . . . , 30 to cos 6'y and i = 31, . . . , 40 to the x distribution. 

The predicted number of events N^^ in the bin i is given by 

iVf = Nb+B{D*'> D\'^)B{D'')TB+Ti , (17) 

where Nb+ is the number of mesons in the data sample and B{D*^ ZJ°7r°) is taken from 
Ref. [T^. In the (JsTtt, e) and (iTvr, /i) sub-samples B{D^) is B{D^ K-71+) [lo|; in {KSn, e) 
and {K3tt,^) B{D^) is Rkstt/k-kB^D^ — > K~7i~^), with Rks-k/Kw a fifth free parameter of the 
fit. Finally, tb+ is the B~^ lifetime [13], and is the width obtained by integrating Eq. [10] 
in the kinematic variable corresponding to i from the lower to the upper bin boundary (the 
other kinematic variables are integrated over their full range). This integration is numerical 
in the case of w and analytic for the other variables. The expected number of events N^^^ 
is related to N^^ as follows 

40 

iV-P = J2 [R^j^jNf) + iVf s'^'^ . (18) 
i=i 

Here, Cj is the probability that an event generated in the bin i is reconstructed and passes 
all analysis cuts and Rij is the detector response matrix, i.e., it gives the probability that 
an event generated in the bin j is observed in the bin i. Both quantities are calculated using 
MC simulation, n]^^^^^ is the number of expected background events, estimated as described 
in Sect. UTTCl 

Next, we calculate the variance of N^^"^ . We consider the following contributions: the 
poissonian uncertainty in Nf"] fluctuations related to the efficiency, which are investigated 
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using a binomial distribution with repetitions and the success probabihty e^; a similar 
contribution based on a multinomial distribution related to Rij ; and the uncertainty in the 
background contribution A^^''^^'^'^. This yields the following expression for af, 



E 



(19) 



bkgrd \ 



The first term is the poissonian uncertainty in A^**^. The second and third terms are the 
binomial respectively multinomial uncertainties related to the finite real data size, where 
^data (-^data) the total uumbcr of decays (the number of reconstructed decays) into the 
final state under consideration [K-n or i^Svr, e or /i) in the real data. The quantities and 
Rij are calculated from a finite signal MC sample (A^mc and A^mc)! the corresponding uncer- 
tainties are estimated by the fourth and fifth terms. Finally, the last term is the background 
contribution o^{N\'^^'^^^^ calculated as the sum of the different background component vari- 
ances. For each background component defined in Sect. IIII CI we estimate its contribution 
by linear error propagation of the scale factor and the error determined by the procedure 
described above. 

In each sub-sample we calculate the off-diagonal elements of the covariance matrix covj^ 
as N'Pij — NpiPj, where Pij is the relative abundance of the bin (^, j) in the 2-dimensional 
histogram obtained by plotting the kinematic variables against each other, pi is the relative 
number of entries in the 1-dimensional distribution, and A^ is the size of the sample. Co- 
variances are calculated for the signal and the different background components, and added 
with appropriate normalizations. 

The covariance matrix is inverted numerically within ROOT and, labelling the four sub- 
samples {Kn or KStc, e or /i) with the index k, the sub-sample functions are calculated, 

Xl = - NnCr/{Nf^ - iVD , (20) 

id 

where A^°^^ is the number of events observed in bin i in the data. We sum these four functions 
and minimize the global with MINUIT (26| . 

We have tested this fit procedure for possible biases using MC samples, generated both 
with and without full detector simulation. All results are consistent with expectations and 
show no indication of bias. 



IV. RESULTS AND SYSTEMATIC UNCERTAINTIES 
A. Results 

After applying all analysis cuts and subtracting backgrounds, 27, 106 ± 367 signal events 
are found in the data. The preliminary result of the fit to these events is shown in Fig. [10] 
and Tab. [ITl The correlation coefficients of the five fit parameters are given in Tab. Illli The 
breakup of the systematic uncertainty is given in Tab. HVl 

As explained earlier, the branching fraction of the decay — > i^'~7r"'"7r"'"7r~ is floated in 



the fit to the full sample. The fit result is in good agreement with the world average [10 
In the sub-sample fits, Rrz-k/k-k is fixed to the value of the full sample fit. 
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1.1 1.15 1.2 1.25 1.3 1.35 1.4 1.45 1.5 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

w, [i] cose , [1] 




-0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

cose^,[i] 



FIG. 10: Result of the fit of the four kinematic variables in the total sample. (The different sub- 
samples are added in this plot.) The points with error bars are continuum subtracted on-resonance 
data. The histograms are the signal and the different background components. The color scheme 
is explained in Fig. [3l 

B. Systematic uncertainties 

To estimate the systematic errors in the results quoted above, we consider uncertainties 
in the following: signal reconstruction efficiency, background estimation, D*^ D^tt^ and 



— > K~7c~^ branching fractions 10|, lifetime [lO(], and the number of mesons in 
the data sample. 

To calculate these systematic uncertainties, we consider 300 pseudo-experiments in which 
13 quantities (corresponding to the above-mentioned contributions) are randomly varied, 
taking into account possible correlations. The entire analysis chain is repeated for every 
pseudo-experiment and new fit values are obtained. One standard deviation of the pseudo- 
experiment fit results for a given parameter is used as the systematic uncertainty in this 
parameter. This toy MC approach also allows the systematic correlation coefficients in 
Table IIIII to be derived in a straightforward way. 

In the following, we describe the parameters varied in the pseudo-experiments: 

• The tracking efficiencies in five bins of slow pion momentum are varied within their 
respective uncertainties. These uncertainties are determined from real data by study- 
ing the decay D*'^tt~^ in the same 140 fb~^ sample as used for the main analysis. 
The tracking uncertainties in different momentum bins are considered fully correlated. 
Therefore, tracking corresponds only to a single parameter in the toy MC, 

• The lepton identification uncertainties are varied within their respective uncertain- 



16 





Ktt P — p 


rfi Ktt f — II 


k'\tv p — p 


P 


i.iyy ± u.i/D ± u.U4o 


i.o(U It u.i/y It u.uo/ 


1 70Q 4- n 1 fi9 4- n nfii 

L.IZO It U.iDZ It U.UDi 


J? M ^ 
Ki{l) 


i.ou/ ± u.ioo ih u.uyu 


i.ODo It U.ioo It U.Uoo 


1 c/in 4- n 071 4- n inn 
i.o4u it u.zii it u.iuy 




n sfia -1- n 00"^ -i- n (T^a 
u.ooo It u.uyo It u.uo^t 


n SQQ 4- n 1 1 n -1- n (v^^ 
u.ooy It u.iiu It u.uoi 


n 4- n 1 Qs -1- n r\A7 
u.ooo It u.iyo It U.U4:/ 


"'iV OTT / I\ 77 


2.072 


2.072 


2.072 




4.91 lb 0.05 lb 0.57 


4.77 lb 0.05 lb 0.56 


4.83 lb 0.07 lb 0.56 




34.3 lb 0.6 lb 2.2 


35.0 lb 0.6 lb 2.3 


36.5 lb 1.0 lb 2.3 


X^/ndf. 


48.3 / 36 


40.6 / 36 


39.6 / 36 




8.3 % 


27.5 % 


31.3 % 





KStt. e = 


Fit to total sample 




1.434 lb 0.209 lb 0.086 


1.376 lb 0.074 lb 0.056 


i?i(l) 


1.813 lb 0.273 lb 0.106 


1.620 lb 0.091 lb 0.092 


^2(1) 


0.764 lb 0.191 lb 0.051 


0.805 lb 0.064 lb 0.036 


RkStt/Ktt 


2.072 


2.072 lb 0.023 




4.83 lb 0.07 lb 0.57 


4.84 lb 0.04 lb 0.56 


Hl)\V,b\ X 103 


34.8 lb 1.0 lb 2.3 


35.0 lb 0.4 lb 2.2 




44.2 / 36 


187.8 / 155 


Px- 


16.3 % 


3.7 % 



TABLE II: The results of the fits to the sub-samples and to the total sample. The first error is 
statistical, the second is the estimated systematic uncertainty. 





-^(1)1^61 


p' 


i?i(l) 


^2(1) 


m\vc\ 


1.000 


0.455/0.399/0.295 


-0.222 /-0.219/-0.179 


-0.054/-0.024/-0.019 


p' 




1.000 


0.648/ 0.413/ 0.540 


-0.889/-0.751/-0.841 


i?i(i) 






1.000 


-0.749/-0.873/-0.763 


^2(1) 








1.000 



TABLE III: The statistical/systematic/total correlation coefficients for jr(l)|V^^| and the three 
form factor parameters in the fit to the full sample. 



ties 



2iM 



The normalizations of the 7 background components are varied within the uncer- 
tainties determined by the background fit in Sect. IIII Ci For continuum, the on- to 
off-resonance luminosity ratio is varied within its uncertainty of 1.5%. Correlations 
between different background component normalizations are taken into account. 

Additionally, we vary the shape in w of the uncorrelated, combinatoric D*^ and fake 
components. We fit each of these distributions by a Lorentz distribution and vary 
each of the parameters within the errors obtained in the fit. 
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fli(l) 


^2(1) 


J'ilWcbl X 103 




Value 


1.376 


1.620 


0.805 


34.98 


4.841 


Statistical Error 


0.074 


0.091 


0.064 


0.37 


0.044 


Tracking 


-0.027 


+0.025 


+0.012 


-1.97 


-0.491 


LeptonID 


+0.012 


+0.024 


-0.011 


-0.39 


-0.096 


Norm - Signal Corr. 


-0.007 


+0.002 


+0.007 


+0.13 


+0.038 


Norm - D** 


+0.005 


-0.023 


+0.002 


-0.04 


-0.041 


Norm - Uncorr 


+0.014 


+0.074 


-0.025 


-0.28 


-0.023 


Norm - Fake £ 


+0.017 


+0.028 


-0.010 


-0.05 


-0.024 


Norm - Comb D*^ 


+0.008 


+0.014 


-0.008 


-0.11 


-0.028 


Norm - Fake 1)° 


-0.009 


-0.014 


+0.007 


+0.06 


+0.020 


Norm - Continuum 


+0.004 


+0.005 


-0.001 


0.00 


-0.003 


Shape - Uncorr 


+0.014 


+0.003 


-0.005 


+0.10 




Shape - Comb 


+0.027 


-0.005 


-0.008 


+0.21 




Shape - Fake 


+0.024 


+0.003 


+0.008 


+0.17 




B{D° Kit) 








-0.32 


-0.089 










-0.82 


-0.227 


B'^ lifetime 








-0.12 


-0.033 


A^(T(45)) 








-0.14 


-0.040 


/+-//00 


+0.003 


+0.006 


-0.003 


-0.15 


-0.043 



TABLE IV: The breakup of the systematic uncertainty in the result of the fit to the full sample. 
The sign + (-) implies whether the fit result moves to larger (smaller) values, if the value of the 
corresponding systematic parameter is increased. 



We vary the fraction /+_//oo = B{T{AS) B+B-)/B{T{AS) B^B^) within its 



uncertainty [10|. This parameter affects the background distributions. 



The uncertainties in B{D*^ D^n^), B{D^ —>■ K^tt^), the number of B~^ in the sample 
and the B~^ lifetime affect only jF(l)|\/c6|; not the form factors. Therefore, they are not 
considered in the toy MC. 

The breakup of the systematic errors quoted in Table [III is given in Table IIVI 



V. MODEL-INDEPENDENT DETERMINATION OF HELICITY AMPLITUDES 

The results discussed in the previous section rely on the form factor parameterization by 
Caprini et al. jsj and it is worthwhile to check this assumption. In this section, we extract 
the form factor shape through a fit to the w vs. cos 9y distribution. For the former quantity, 
we consider six bins in the range from 1 to 1.5. For the latter, we have also six bins between 
— 1 and 1. The contribution from events with w > 1.5 is fixed to the very small values 
predicted by the results of the parametrized fit. 
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A. Fit procedure 



From Eq. [10] we can obtain the double differential decay width d'^T / dw d cos 6v by inte 
gration over cos and x- 
If we define 



(21) 



and 



1 — 2wr — r I w — 1 



w — 1 



7°°(u;) = V^I^iw + l)'hl^{w)\V,,\' 1 + ^ (1 - i?2H) \ , (22) 



2 



the double differential width becomes 
dwd(cos6'y) 3 



Fp (sin^ Ov (7++ + 7—) + 2 cos^ Oy 7°°) ■ (23) 



The quantities 7^^*^ correspond to the w-dependence of the different helicity combinations, 
times kinematic factors. The one dimensional distribution, as given in Eq. [TT], depends only 
on the sum of these three combinations, 

^ d^;^ = — (7++ + 7-' + 7°°) • (24) 

The bin contents of the two dimensional histogram in w vs. cos^y can be obtained by 
integration of Eq. [23] over the corresponding bin area and considering the reconstruction 
efficiencies and detector response as described in Eq. [18] Each bin content can be given as 
the linear combination of two linearly independent parts. The integration of the angular 
distributions is executed analytically, as described in Sect. lIIIE| the integration with respect 
to w defines a set of free parameters, 

rf = [ dw7^^ (25) 



where Wj = {wi, W2, ■ ■ ■ , wj} = {1, 13/12, . . . , 1.5} are the bin boundaries of the 6 bins in w. 
Additionally we define 7^ = 7"^+ + 7 and Tj = T^^ + T'- . 

We define the function x^,^^ = Ei=i Ej=i [ ^'aj^e.p J > ^^ich depends only on the 

parameters Tj and r°°. Here A^°'''' gives the number of events observed in on- resonance 
data, N^^'P the number of expected events, as defined in Eq. [18], and a^exp the uncertainty in 
the expected number of events, as given in Eq. [19] This expression is minimized numerically 
using MINUIT to determine the partial integrals. 

Investigation of dedicated sets of toy Monte Carlo showed the reliability of the procedure 
described above. The results show no indication of bias in either the mean or the errors. 
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FIG. 11: Results of the fit of the helicity amphtudes (red crosses) compared to the prediction 
obtained by using the parametrization prescription by Caprini et al. 0] (sohd black line). The left 
plot shows the results for T^^, the right one for Vf. Only the statistical error is shown. 

B. Results 

Given the very high amount of background in the K~7i^7i~7i^ mode, we use only 

the Ktt, e and Kn, fi channels to determine the partial decay widths for each of the helicity 
components. Tables IVl and I VII give the results of the fits, where the systematic errors quoted 
in these tables stem from the same sources as given in the breakdown in Table IIVI It is 
dominated by the track reconstruction errors of the three charged tracks i, it, K and the 
uncertainty of the tt^ reconstruction. The of the fit is in good statistical agreement 
with the number of degrees of freedom, we obtain x^/ndf = 82.6/60 or a probability of 
P^2 = 2.8%. The results are shown in Fig. [TTl There is good agreement with the results of 
the parametrized fit. 

VI. SUMMARY AND DISCUSSION 

We have reconstructed about 27,000 S+ D*°i+Ui decays in the 140 fb"^ of Belle 
T{AS) data. A fit to the theoretical expression for the four- dimensional differential decay 
width (Eq. [TO!) . assuming the parameterization of the helicity amplitude given by Caprini 
et al. yields a measurement of \Vcb\ times the form factor normalization at zero recoil, 
^(1)1^61 = (35.0 ± 0.4 ± 2.2) X 10"'^. At the same time we determine the parameters of the 
Caprini et al. parameterization, = 1.376 ± 0.074 ± 0.056, i?i(l) = 1.620 ± 0.091 ± 0.092, 
R2I1) = 0.805 ± 0.064 ± 0.036. The branching fraction of the decay B+ D*°i+ue is 
measured to be (4.84 ± 0.04 ± 0.56)%. For all numbers quoted here, the first error is 
the statistical and the second is the systematic uncertainty. All results are preliminary. 
These measurements are in agreement with previous investigations of the decay — >■ 
D*'i+ue 0,0,8- 

A direct, model- independent determination of the form factor shapes has also been carried 
out and shows good agreement with the HQET based form factor parametrization by Caprini 
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TABLE V: Obtained results for , compared to the central values of the parametrized fit. 
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( 2.421 ± 0.175 ± 0.258 )xlO- 


-4 


2.483 xlO-^ 




pOO^ 


w 


G(i|,1.5) 


( 2.683 ± 0.186 ± 0.298 )xlO- 


-4 


2.514 xlO-^ 





TABLE VL Obtained results for F™, compared to the central values of the parametrized fit. 
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